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Abstract 

The 1-loop quantum weight of the SU(N) KvBLL caloron with nontrivial holon- 
omy is calculated. The latter is the most general self-dual solution with unit topo- 
logical charge in the 4d Yang-Mills theory with one compactified dimension (finite 
temperature). 

1 Introduction 

^ ■ The finite temperature field theory is defined by considering the Euclidean space-time 
which is compactified in the 'time' direction whose inverse circumference is a tem- 
perature T, with the usual periodic boundary conditions for boson fields and anti- 
periodic conditions for fermion fields. In particular, the gauge field is periodic in time, 
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so the theory is no longer invariant under arbitrary gauge transformations. Only time- 
periodic gauge transformations are allowed and hence the number of gauge invariants 
increases. The new invariant is the holonomy or the eigenvalues of the Polyakov line 
that winds along the compact 'time' direction [1]: 



X 



L = P exp 




(1) 



This invariant together with the topological charge and the magnetic charge can be 
used for the classification of the field configurations [2] , its zero vacuum average is 
one of the common criteria of confinement. 

A generalization of the usual Belavin-Polyakov-Schwartz-Tyupkin (BPST) instantons [3] 
for arbitrary temperatures and holonomies is the Kraan-van Baal-Lee-Lu (KvBLL) 
caloron with non-trivial holonomy [4, 5, 6]. It is a self-dual electrically neutral con- 
figuration with unit topological charge and arbitrary holonomy. This solution was 
constructed by Kraan and van Baal [4] and Lee and Lu [6] for the SU(2) gauge group 
and in [5] for the general SU(N) case; it has been named the KvBLL caloron (recently 
the exact solutions of higher topological charge were constructed and discussed [7, 8]). 
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Figure 1: The action density of the SU(3) KvBLL caloron as function of z,t at fixed 
x = y = 0, eigenvalues of A4 at spatial infinity are ]i\ = — 0.307T, ^2 = — 0.013T, ^3 = 
0.32T. It is periodic in t direction. At large dyon separation the density becomes static 
(left, Qi f 2 = 1/T, £3 = 2/T). As the separation decreases the action density becomes 
more like a Ad lump (right, Qi r 2 = 1/ (3T), ^3 = 2/ (3T)). The axes are in units of 
inverse temperature 1/T. 

There is a plenty of lattice studies supporting the presence of these solution [9], see also 
[10] for a very brief review. In a recent paper [11] the caloron ensemble was studied 
analytically, although some contributions were neglected there, the results are in very 
good agreement with phenomenology. 

The holonomy is called 'trivial' if the Polyakov loop (1) acquires values belonging 
to the group center Z(N). For this case the KvBLL caloron reduces to the periodic 
Harrington-Shepard [12] caloron known before. The latter is purely an SU(2) configu- 
ration and its quantum weight was studied in detail by Gross, Pisarski and Yaffe [2]. 
The KvBLL caloron in the theory with SU(N) gauge group on the space R 3 x S 1 can 
be interpreted as a composite of N distinct fundamental monopoles (dyons) [13][14] 
(see fig. 1 ). As was proven in [5, 17], the exact KvBLL gauge field reduces to a su- 
perposition of BPS dyons, when the separation Qj between dyons is large (in units of 
inverse temperature). On the contrary, the KvBLL caloron reduces to the usual BPST 
instanton, when the distances Qi between all the dyons become small compared to the 
inverse non-triviality of holonomy. 

We refer the reader to the papers [5, 17] for the detailed discussion and construction of 
the caloron solutions, to the original works [4] for the SU(2) case and to further works 
on higher topological charge solutions [7, 8]. 

This paper is in the series of papers [15, 16, 17, 18, 19] where we calculate the func- 
tional determinant for KvBLL calorons with nontrivial holonomy [4, 6] in the finite- 
temperature Yang-Mills theory. 

Here we calculate the 1-loop gluonic and ghost functional determinants for the case 
of an arbitrary SU(N) gauge group. The calculation is performed in the limit of far 
separated dyon constituents and up to an overall numerical constant. The constant for 
the gluonic determinant remains known only for the SU(2) case [15] 2 . 
We find new 3-particle interactions arising between constituent dyons due to the glu- 

2 In the previous paper [17] we have proved that the corresponding constant is zero for the 
fundamental-representation determinant, but for the adjoint representation the constant is not expected 
to vanish. 
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onic determinant. These terms were not present in the fermionic (fundamental repre- 
sentation) determinant and also vanished in the SU(2) gluonic determianant. 
Similar to our previous results, the determinant is infrared divergent, the leading di- 
vergence is proportional to the volume or the system, and there are of course several 
subleading divergent terms. It is not surprising and was known long ago [25, 2] that 
non-trivial holonomy increases the effective action by a factor proportional to the vol- 
ume. Nevertheless that does not make the studies of non-trivial holonomy unphysical, 
since in the ensemble of many calorons the moduli space integrals can compensate the 
above divergences. 

Since there are subleading divergences coming from the Coulomb tail of dyon fields, 
it is natural, that the result would also depend on the position of the large ball, with 
which we make the infra-red cutoff. We will display this dependence, but we note that 
it is unphysical unless the box in not a real border of plasma region. One could also 
expect that the quantum corrections will dump the Coulomb tails of dyons. 
We present the relevant notations in section 2 and illustrate the notations by the old 
results. The method of computation is described in 3 and the actual computations are 
carried out in subsequent sections and appendices. The final result is presented in 
section 6. 

We do not draw here any physical conclusions on the behaviour of the whole caloron 
ensemble since that is now a separate business [26], [11]. Our results could be useful 
for extending the works [26] to the SU(N > 2) case and including the corrections due 
to the non-zero modes to the work [11]. 

2 Notations and Review 

Consider the SU(N) YM theory and a caloron solution with the asymptotics 3 

A H ^>°° InSjA diag(/*i,...,/i N ). 

For the SU(2) case the standard choice is }i\ = —co = co where < to < \. 
As usual, we set the temperature T = 1 throughout the computation, and restore the 
temperature dependence only in the final result. 

The caloron can be viewed as composed of dyons (BPS monopoles with A4 playing the 
role of a Higgs field), the inverse dyon size v being defined as 

v\ = m +1 - m ; v N = ]i x - }i N + 1. (2) 

Traditionally the first N — 1 dyons are called the 'M - dyons' and the N~ th dyon is 
called the 'L - dyon', because an additional gauge transformation is need for it to have 
correct asymptotics. 
We also introduce a notation 

v mn = 2n{]i m - ]i n ) mod 2n. (3) 
which coincides with v = V21 and v = n\i used previously in the SU(2) calculations. 

3 We use notations consistent with [5]. 
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The positions of dyon centers are denoted by y ; . The distance from the i th dyon center 
to a point x is denoted as x — 



ra for the SU(2) 



case the standard notation is 



r i = s ; ri = r [15]. The distance between dyon cores is denoted by r 
It is convenient to use a so-called 'algebraic gauge', in which the asymptotic gauge 
field is vanishing at the expense of introducing twisted boundary conditions for field 
fluctuations. The twist a(x,l/T) = e~ lT a(x,0) is hence related to the holonomy as 
t = 2;rdiag(^i, The holonomy and, correspondingly, the twist could also be 

multiplied by elements of the center of the gauge group e 2/n w. It does not affect the 
adjoint gauge field and determinant but it affects fundamental determinants [17]: 

logdet(-V^) = E (>»( T „ )r „,„_ 1 + Ip (T „)^)-^-fcl 



+c N + -logp + 0(l/r) 



(4) 



where 



cn 



13 TT 2 log 7T £'(2) 



72 216 



+ 



P is a periodical function with a period 2n such that 



(5) 



N q 2 (2n-q) 2 
P(v) = ^ 127r2 ; <7 = vnod27r 



(6) 



Determinant in the adjoint representation of SU(2) reads [15, 19] 



logDet(-Djj) = VP(v) + 2nP"(v)r l2 + 

2 5 
+ 3 lo gM+ 3 lo g( 27r ) + c 2 



37T-4V 
37T 



logv + 



+- 



1 1 23/r 87 £ 74 4 , 
+ = + ^-^-^-^log 



v v 54 



3n 9n 3n 



w r 



12 



71- 



3tt-4v 
3tt 



logv 



+ o 



12 



(7) 



Now we proceed to the calculation of SU(N) adjoint-representation determinant. 



3 Method of computation 

For self-dual fields the gluonic and ghost determinants over non-zero modes for the 
background gauge fixing are related [27] to the adjoint scalar determinant in the same 
background: Det'(W„ v ) = Det(— D 2 ) 4 , where W„ v is the quadratic form for spin-1, 
adjoint representation quantum fluctuations and D 2 is the covariant Laplace operator 
for spin-0, adjoint representation ghost fields. So the total contribution to the effective 
action of gluon and ghost determinants is 21ogDet( — D 2 ) which corresponds to two 
physical degrees of freedom. 

We calculate the quantum determinant by integrating its variation with respect to pa- 
rameters V of the solution, following [17, 16, 15, 28]. In this case the problem reduces to 
four dimentional integral of the gauge field variation multiplied by a vacuum current, 
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which can be expressed through Green function known implicitly for any self-dual 
configuration 

8 logDet(-D*[A]) = _| ATr(w) (8) 

where / ? , is the vacuum current in the external background, determined by the Green 
function: 

] ll =B ]l Q + QD ]l . (9) 

Here Q is the periodical Green function of the covariant Laplas operator in adjoint 
representation 



-D 2 x G{x,y) = 6®{x-y) (10) 
+00 

S(x,y) = E G(x 4/ x;i/4 + n / y). (11) 



The Green functions in the self-dual backgrounds are known explicitly [29, 30] if the 
gauge field is expressed in terms of the Atiyah-Drinfeld-Hitchin-Manin (ADHM) con- 
struction [31]: A ?/ = v f d^v. These look quite simple for the fundamental representation 
[32] 



G {X ' y) ~ ±n 2 (x-y) 2 ' (12) 



but become more complicated for the adjoint representation [33, 34, 35]: 

G a \x,y) 



\Trt«{v{x)\v{y))t h (v{y)\v{x)) 



4n 2 (x-y) 2 
1/2 

+ - — j I dz\dz2dz^dz^M{zi,Z2,zo ) ,z^) 



(13) 



-1/2 



x lTr(v + (x,z 1 )V(x,z 2 )^)Tr(v + (y,z 4 )V(y,z 3 )^) , 



where t a are Hermitian fundamental-representation generators of SU(N) normalized 
to \xt a t b = \5 ab ; V(x,z) is one of the components of v (see eq.(65)) and M is a piece- 
wise rational function 4 . Fortunately we do not need an explicit form of this function 
for the SU(N) caloron since in the large separation limit the contribution of the last 
term (or "M-term") is exponentially small away from the dyons. Near the dyons the 
field is essentially reduced to the SU(2), so one can use there the results of [34, 15]. 
In what follows it will be convenient to split the periodic propagator into three parts 
and consider them separately: 

Q{x,y) = £ r (x,y)+£ s (x,y)+£ m (x,y), 

g s (x,y) = G(x,y), (14) 
Q T {x,y) +Q m (x,y) = G(x 4 ,x;y 4 + n,y) , 



'see Appendix B to [15] for its explicit form in case of SU(2) 
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here Q m (x,y) coresponds to the part of the propagator arising from the M-term. The 
vacuum current (9) will be also split into three parts, "singular" , "regular" and "M", 
in accordance with (14) 

h = Jl+ j; + j?- (is) 

As was proposed in [15] we divide the space into regions surrounding the dyons and 
the remaining space (outer region). Near each of the dyons the gauge field becomes es- 
sentially the SU(2) dyon configuration plus an additional constant-field background. 
In this region we can use the results of [15]. In the outer region, far from the exponen- 
tial cores of the dyons, the vacuum current considerably simplifies and we only have 
to perform integration in (8). 

In the following two sections we give results for these two domains and in section 5.1 
we combine them together and integrate over the space. 



4 Core domain 

In this section we write a contribution to the variation of total determinant arising 
from the core region of a dyon. We take a ball of radius R around the dyon. In that 
region the field is approximately the one of the Slf(2)-dyon, embedded along one of 
the simple roots, plus an extra constant A4 field [5]. More precisely in the fundamental 
representation the gauge field near the Z th dyon is a zero N x N matrix with only 2x2 

block at r*h position filled by the BPS dyon gauge field, plus a constant diagonal N xN 
matrix [17] 

< block 2x2 = A dyo >Z/ _ ^ + ^ (^1+^ s ^ hx2 f (16) 
A out S1 de & block 2x2 = ^ • diag ^ _ _ ^ g ^ 

Under the action of the SU(2) sub-group, the adjoint representation of SU(N) splits 
into one triplet, 2(N — 2) doublets and (N — 2) 2 singlets. The determinant of arbi- 
trary SU(2) configuration embedded into SU(N) is then expressed as a sum of SU(2) 
adjoint-representation determinant plus 2(N — 2) fundamental-representation deter- 
minants 5 . 

As is seen from eq.(16), the SIT (2) dyon field is accomplished by the constant diagonal 
matrix. This matrix can be killed by a gauge transformation, which is not periodical, 
and thus can change the determinant. It is equivalent to the additional twist of bound- 
ary conditions for the 2(N — 2) fundamental representation determinants. 
As a demonstration let us consider the SU(3) case. The fundamental gauge field reads 



A ¥ = \ 


( A dyon 
A fi 2x2 








V 







+ 2mSu4\ ™ (17) 



'See for example [36] where that was done for the instanton solution. 
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in the adjoint representation of SU(3) in an appropriate basis it becomes 

0\ 



1 3x3 



A, 







V o 







-A 



dyon 
\i 2x2 







A 



dyon 
ft 2x2 



(18) 



/ 



So that there is one block 3x3 giving adjoint representation dyon field and two 2x2 
blocks, giving fundamental representation dyon accomplished by a unit matrix. As it 
was shown in [17] this extra unit matrix changes only an IR divergent part of the dyon 
determinant (the one, depending on radius of the ball). These divergences cancel with 
the terms in the outer-region determinant depending on the radius of the holes R, as it 
is shown in Appendix A). So we can freely drop them. 

Summing up one adjoint (42) and 2(N — 2) fundamental dyon determinants (41) for 
all N dyons we obtain the following contribution to the derivative of the Caloron de- 
terminant from the considered domain 



(6 + N)v n log(v„) 



+ IR, 



(19) 



where "IR" denotes the IR divergent terms. 



5 Outer domain 

We proceed to consider the far domain, i.e. the region of space outside dyons' cores. 
Caloron field becomes diagonal with 0(e~ v ' r ') precision and this simplifies signifi- 
cantly the results. For instance the 4 component of the fundamental caloron gauge 
field reads 

Ar = /r»(2^ + A--^J_) (20 ) 

In what follows we will consider the derivative of the determinant with respect to }i m . 
It turns out that in this domain only A4 depends on ]i m nontrivially [5] . Thus we need 
only 4 th component of the vacuum current as it follows from (8). As we know from the 
SU(2) case this component of the current is especially simple [15] 

JT {2) = ^T3P'(v+---) (21) 
the natural generalization of this expression is 



t ^ N 

/f (n) = ^ag 

n,m=l 



. F > 2n( Fm - Fn ) + — -- — + 



2 \ 2v m 2v m —\ 2v n 2t n . 



(22) 



The expression in the brackets is simply the eigenvalue of the the gauge field (20) in 
the adjoint representation. This formula is definitely right for large r m , where the field 
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becomes almost constant [2, 37] and generalizes the SU(2) expression. Moreover we 
check it by a direct computation in Appendix B. We conclude that 



-fxidpAM = \E d rP (v mn + ^r- ^— - + (23) 

z n,m V Lr m Lr n ^'n-lj 

where 

v mn = 2n(}i m - ji n ) (24) 

The variation over V can be integrated up to a constant, and the integral over space 
would be performed in the next section. 



5.1 Integration 

In order to get a variation of the determinant we have to integrate in eq.(23) over the 
space with N spherical holes of radius R . The following integrals will be very helpful 

1 1 1 1 \ 2 j3 

1 d°x ~ 

2r m 2r m -i 2r n 2r„_i/ 

Tt {jm,n ~l~ Vm—\,n—\ ^m,n—\ ^m—\,n "4" ^m,m—\ ^n,n—l) (25) 



1111 



+ — - — ) d 3 x 



2t m 2f m —i 2t n 2t n —\ 
-3/rlog — : ^ '- '- (26) 

Vn,m-1 Tn,m,n—\ ^m,n—\,m—\J 

where r nm = \y n — }/m \ is a distance between dyons and 

2 Y\mn = r lm + r mn + r nl (27) 

is the perimeter of the triangle, formed by Z -th , m _th , n _th dyons. Sign ~ means that we 
drop all the terms dependent on the radius of the holes R since they cancels precisely 
with dyons IR divergences as discussed in Appendix A. To derive the last equation we 
used 

/ d 3 x ~ — An log r nm i + C. 

J T n T m Ti 

It is important to point out that the eq.(26) is not applicable for the case m — n ± 1, 
since it diverges. The reason is that the divergences near dyon cores are not balanced 
anymore. Nevertheless it is straightforward to verify that if one replaces a zero r nn 
under logarithm in eq.(26) by some fixed e, then it is is still valid up to a constant, 
which cancels in the final result. 
So we can integrate in eq.(23) 

1 E^f v - + ^-^ J --^ + ^ J -V 3 ^ ( 28 ) 

m,n 

n 



2 m;n V 2f m 2r m _\ 2f n 2r w _i, 

/ , ~7~P ( v mn) (?m,n "I" ?m—l,n—l ~ ?m,n—l ~ ?m—l,n r m,m—l "I" ?n,n—lj 



4 

m,n 



v— i [ v mn\ TC , f r m,n—l r m,n,m—l ^n,m—l,n—l\ . V -1 ^ n/ \t^(3) 

L 27r lo § ( r ir ir i i J + ^ 2 P(Vww)y 

m,n \'n,m— 1 'n,m,n— 1 ' m,n— \,m— 1 J n,m 



We denote [v mn ] = v mn mod 2n. To simplify the above expression we use the identity: 

V -1 \ymn\ 7T i ^m,n,m—\^n,m—\,n—\ v -1 i Vi /<->r>\ 

L — 17 — lo § 7"^ — 7^ — = L Vn lo § r «,«,«-i - L lo S r n,n-i- (29) 

m,n 'n,m,n—Vm,n—\,m—\ m ,n n 

Then eq.(28) becomes 

log det(- D ^ = / 1 £ P (v„„ + ±. _ _L_ _ _L + A £ 



(30) 



l v m«l / T m n —\ \ r— i . r— i 

w ~ log \ ~l^2v n log r m ,„, m _i + 2 log r W/M _! + 

m,« 171 Vn,m-1/ m/W M 

y~! T^ 5 ( v mn) {jm,n + r m—\,n—\ ~ r m,n—l ~ r m—l,n + r m,m—l + r n,n—l) + 

+ Ei p ( v -)^ (3) 

The "R-terms" are exactly the ones of eq.(44) but with R standing as a lower limit of 
integration, this provides the cancellation of them when we add the core contribution. 
The second equality in (30) is valid when the variation does not involve changing of 
the far region. Note that the correction comes only from the far region. So we can 
calculate it. It comes from the next P"" term in the Taylor series, this term obviously 
involves 4-center Coulomb integrals: 

^ (3D 



rir 2 r 3 r 4 



taken over R 3 with holes around centers. Since this integral converges both in the IR 
and UV (near the holes) , it can involve only logarithms of some dimensionless combi- 
nations of distances between these four points, divided by the distance. In the approx- 
imation that the dyons are spread homogeneously, such terms would be of order of 
unity and we neglect them. The only large logarithms come from the case where three 
of the four points coincide, in this case the integral diverges as logarithm near the z-th 
dyon: 

1 log(r„/R) 

4tt ) + P(l/r f/ ) (32) 



R4\B R 



i'i 



So for the correction to logDet(— D 2 ) one sums all the contributions of the form (32). 
Note that P IV = is a constant, so some terms cancel in the sum. The result for N > 2 
is: 

logDet(-D ) correction = 7 L ' ( 33 ) 

b7r n>m r nm 

for N = 2 the coefficient is doubled and becomes — j^, since there are more coincident 
points. It matches our SU(2) result (eq. (60) in [15] ). 
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6 The result 



From eqs. (19,30) we can conclude that for large dyons' separations, q m <C l/v m + 
l/v m _i, the SU(N) caloron determinant is the sum of these expressions plus some in- 
tegration constant and improvement (33). Restoring the temperature dependence 
we obtain 



logdet(-D^) 



(34) 



m,n 



[V mn \ 7t I T m>n _i 

2n ° g I r n , m -i 



E 2v " lo 8 

'm,n,m—\ ' w,n— 1 



m,n 



+ E ( Vmn )T ( r m,n + r m—l,n—l r m,n—l r m—l,n + r m,m—l + r n,n—l) 

m,n 

+ E ^(v„,„)tVct - e (6+w) :"' ogl/ " + Eiogv„ - ^-^f + 



m,n 



c N 



Note that the coefficient — should be doubled for N = 2 case. 

The contribution to the effective action from non-zero modes of gluons and ghosts 

would be 

Det(-D 2 ) , „ . 2 



SS 



eff 



log 



(Det'W„v) 1/2 



logDet(-D 2 



(35) 



The constant will of course contain a standard UV-divergence cjv = c + y log }ipv, 
coming from the instanton determiant [38], where }ipy is a Pauli-Villars mass. This di- 

4N 



vergence, together with 



coming from zero modes, gives the standard Yang- 



Mills /3-function and is commonly incorporated into the running coupling: 



jipy e 



T" <? S 2 (}i PV ) 



(36) 



where A is the scale parameter obtained here through the 'transmutation of dimen- 
sions'. 

Now let us combine with the result for SU(N) caloron zero modes [23, 22] and the 
classical action 8n 2 /g 2 (}ipv)- The caloron measure is [23] 



W ^ 2 6N 7T 4N 



! + Ei — — + — 



(37) 



So the total contribution of one caloron to the effective action becomes 



gSeff 



8n 2 



2N 



X 



X 



i + E 



1 1 

+ 



(38) 



We have collected the factor / kne~ ryE T/ A because it is the natural argument of the run- 
ning coupling constant at nonzero temperatures [39, 37]. When we have done so in the 
SU(2) case [15], we have got a constant numerically very close to 1, so we expect the 
constant C to be of order of unity. 
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A Cancellation of IR divergenses of dyons 

Consider the field near the dyon constituent of SU(N) caloron. In the fundamental 
representation it is given by the eq.(16). In the adjoint representation this field looks 
like one block with the SU(2) BPS dyon in the SU (2) -adjoint representation (3 x 3) and 
2(N — 2) blocks with the SU(2) BPS dyon in the SU (2) -fundamental representation 
(2 x 2) plus a constant part, specified below: 

-iA a J nst = 7rdiag(0 / / / -(/ii + /*2-2/i 3 ) / -(/ii + /i2-2^3), 
+ (..), + (..), + (ui + }i2 - 2u N ), +(u 1 + u 2 - 2u N ), 
,/^3 — f*4,^4 — ^3, ...(all pairs without ui,U2)...,0,0, ...) (39) 

There are N — 2 zeroes in the end, corresponding to elements of the Cartan subalgebra 
other than T 3 . Let us check the size of the matrix (39): 



3 + 4(N - 2) + (N - 2)(N - 3) + (N - 2) = N 2 - 1 
as it should be for the adjoint representation of SU(N). 

This constant background is exactly equivalent to twisting the boundary conditions. 
So we have to sum logarithms of determinants for one adjoint-representation SU(2) 
dyon, 2(N — 2) differently twisted fundamental-representation SU(2) dyons and a de- 
terminant for (N — 2)(N — 3) different constant A4 field eigenvalues. It would be in- 
teresting to check that this is asymptotically the same as the adjoint determinant for 
the far region that we will calculate below. 

In [17] we proved that the twisting of boundary conditions for the fundamental repre- 
sentation determinant results in shifting of the argument of P, where P is the standard 
perturbative potential 



(40) 

mod In 



Now we can just write the result (for the first dyon, for simplicity), as a sum of the 
formula taken from [15] for the triplet (SU(2) adjoint representation): 

9plogDet(— D ) 

near dyon 

= d v I c dyon v x - |vi log(vi) + logv m + J P (^2nv 1 - Anr 2 dr J (41) 
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with 2(N — 2) formulas from [17] for the SU(2) dyon with twisted boundary condi- 
tions ("twist" is a corresponding matrix element of eq.(39)) 6 : 

2(N-2) , \o (y \ 

d v logDet(- V 2 ) near dyon = E d T |%o« y i - ° g 6 1 y \ 
+ J \ [p (27t( H i - (Mi + p 2 )/2) + ^ - 



+P [ 2n{ n - (hi + m)/2) - Tr ~ i{A a c d J nst )i l+ 2 



4nr 2 dr 



(42) 



and with a constant-field determinant for the remaining matrix elements of — D 2 ( i.e. 
twists without background field) 



(N-2)(N-3) R 

aplogDet(-D 2 ) n - st dy0 n = E ^ / 2 ? (" Z ' (A ^ )2+4 ( N - 2 ) + V Wdr - 







Totally we get: 

3p logDet(-D 2 ) near dyor 
+ terms 1 



dp I c dyon v n 



8 N-2 
3 + — 



log(i/ m )v m +logv„ 



(43) 



And it is easy to check explicitly that the "R-terms" exactly match the asymptotics of 
far-from-dyons domain (see the Section 5) 



(R-terms) = lim / Y SP 



1 1 

in{m - m) + — 



i 



2r t 2r l _ x \2r ] 2r y _ a 



47rr 2 dri 
(44) 

So we conclude that the "R - terms" are trivial and exactly match the contributions 
from the outer region, as it should be, of course, since the result cannot depend on the 
radius of the auxiliary balls that we have chosen. 



B Calculation of the currents for outer domain 

B.l Singular current 

The contribution of the singular part of the propagator to the variation of the determi- 
nant is 4 times the fundamental representation result [27, 2] , if we write this variation 
in terms of fundamental representation. So we just take our old result from [17] (that 
formula was not written there explicitly, it was in our intermediate computations). For 
the component } sl it is quite natural to introduce bipolar spatial coordinates with unit 

6 Note that to get the normal SU(2) dyon we have to center the interval (^1,^2) at zero, because in 
SU(2) we obviously have the zero-trace condition \i\ + ]ii = 
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repers ?, ■ = ^ , §j = , n^, = -^ffy. In these coordinates the current (already multi- 
plied by 4) is 

n = (45) 



z (n + s f ) y/- (gj - Tj - Si) (d + rj - Sf) (gj - rj + sjj (gj + rj + sjj 
v 47r 2 r 2 s 2 ( e + n + SiY 

1% = ; J% = (47) 

We also remind the notations: Y\ = x — i/j is a vector from the z'-th dyon center to the 
current point, s ; = r,_i; and Qj = \yi — T/j-i| is a distance between these two dyons. 
Also standard "circular rule" r^+i = ^i is implied. 

B.2 M-term current 

Let us prove that the contribution to the current from the M-term of adjoint propagator 
is zero with exponential precision (i.e. it decays exponentially out of the dyon cores). 
As was shown in [15], when making the propagator periodic the M-term simplifies to 

1/2 

Q mab (x,y) = / dzdz'M{z,z')Tr (v 2 \x,z)v 2 {x,z)T a ^ Tr (z; 2+ (y,z> 2 (y,zV) , 

-1/2 

(48) 

since the property v{y n ,z) = e 2mnz v(y,z) used to derive that result, still holds for the 
SU(N) ADHMN construction. Here M(zi,z 2 ,z,z) = 5{z\ -z 2 )M(zi,z). 
First of all we note that only the lower components of v are left and only the Cartan 
(diagonal) components are nonzero: From eq.(78) we see that for each m the function 
s m (z) and hence v 2 n (z) is peaked near z = }i m and exponentially decays away from 
this point. So v 2 ^(x,z)v 2 (x,z) ~ Smn with exponential precision. This leads us to 
conclusion that 

g mab (x,y) oc s aeCartan S beCartan , g mab {x,y) = g mab (y,x) . (49) 

The second equation means that the terms with derivatives in the expression for the 
current (9) cancel each other. It follows from the first one that the adjoint action of A 
on g m gives zero since both lie approximately in the Cartan subalgebra. Therefore we 
conclude that 

/ f ? - . (50) 

B.3 Regular current 

The adjoint-representation regular current is 

J ab = D a x c g cb (x,y)+g a c(x,y)Df (51) 
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where a, b, c = 1..N 2 — 1 and we take the regular part of the propagator: 



m ab (*,y) = E 



^ Q 8n 2 (x-y n ) 2 



Tr 



t«{v(x)\v(y n ))t b (v(y n )\v(x)) 



y l n 



(52) 



It is possible to rewrite these formulae in the fundamental notations and evaluate ex- 
plicitly. Some details of the calculation together with a short review of ADHM con- 
struction are presented below. We denote the adjoint indices by a, b, c = 1..N 2 — 1 and 
fundamental indices by k, l,m,n = 1..N. 

First we represent the covariant derivatives in the fundamental representation. With 
the help of the identities 



tr \t At B 
tr ( t a At d B ) D 



tr 
tr 



t a (D ¥ A)t h B-t a At h (BD ¥ ) 
t a (ADu)t b B-t a At b {DuB) 



(53) 



2tr\t a v v t v v 1 

one gets for J the obvious four terms plus 5^ ^ 3 y from the derivative acting 

on the denominator. 

All the terms in the adjoint current have the form tr[t a Bt C]. The variation of the de- 
terminant has the form —5A c T c ah J ab . To write it in the fundamental representation we 
use the identities: 

where A = AH 1 and 



AX b = 2tr(t b [t\A}) 
^ = l/2(S tt S jk -l/NS ij S i 



(54) 



kh 



(55) 



The Hermitian generators t a are normalized as tr(t a t b ) = 1/2 (for SL7(2) these are 
tV2). We get: 



SA c T c ab tr(t a Bt b C) = -tr(B) tr(M C) - hr(C) tr(M B) 



(56) 



So in terms of the fundamental indices = 1...N) we get for the current J l J (that is to 
be coupled to A in the fundamental representation to get the variation of determinant) 



Here 



E 



n^0,{B,C} A7l2 ( X V n ' 

-Si* E 



tr(B) & - tr(C) B 1 ' 



tr(E) pi - tr(F) E*> 



vlVy n 



b ; F = vt v x 



(57) 



(58) 



and we put y = x (so that now y l n = x 1 — 5 lA n ) according to eq.(9). The set {(B, C)} 
consists of 4 pairs taken from eq.(53) and eq.(51): 

{(B,C)} = {(D x vlvy nl vlv x ), (vlv Vn ,D x vl n v x ), {v\v yn D y ,v\v x ), {v\vy nl v\v x Dy)) 

(59) 
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Since the current and the field are approximately diagonal in the fundamental repre- 
sentation, we consider the diagonal components (J^) 1 = (Ji) tt - In these notations the 
contribution to equation eq.(57) can be rewritten as 

(/*,)'(*) = E E ( E fa,/* y (B, Q - q B.) - ^ (£; Fj - F, E.) J 

(60) 

From eq.(59) and eq.(60) we get 8+2 terms in the resulting contribution to the current. 
Now in order to calculate explicitly these B, C, E, F we need ADHMN construction. A 
brief review and a calculation follows. 



B.4 Expressions of the ADHM construction 

The basic object in the ADHMN construction [31, 30] is the (2 + N) x 2 matrix A lin- 
ear in the space-time variable x and depending on an additional compact variable z 
belonging to the unit circle: 

K j Aj^z) , K = m, 1 < m < N, 

V z ' x ) = I (B( Z )- W )J , K = N + a, l<a<2, (61) 

where a, ft = 1,2 and m = 1, . . . , N; cr„ = (ur, I2). As usual, the superscripts number 
rows of a matrix and the subscripts number columns. The functions A? (z) forming 
a N x 2 matrix carry information about color orientations of the constituent dyons, 
encoded in the N two-spinors £: 

^{z) = 5{z- } i m )q. (62) 

The quantities £g transform as contravariant spinors of the gauge group SU(N) but 
as covariant spinors of the spatial SU(2) group. The 2x2 matrix B is a differential 
operator in z and depends on the positions of the dyons in the 3d space y m and the 
overall position in time £4 = X4: 

6%d z A%(z) 

with 

A(z) = A F cr F , A(z) = 2niy m {z), A 4 = 2m £4, (64) 

where inside the interval ji m < z < }i m +i, y(z) = y m is the position of the dyon 

with inverse size v m = }i m +i — }t m - 

One has to find N quantities v„(x), n = 1...N, 

v K (x)-t V « 1 ^ ' K = m ' 1 < m < N ' ((ft 
V »W ~ \ v*(z,x) , K = N + oc, l<a<2, 

which are normalized independent solutions of the differential equation 

1/2 

X f l(z)vl m + [B\z) - xrf$ff{z,x) = 0, vj lm vl l + J dzvl 2m v 2 n a = 8%, (66) 

-1/2 



15 



or, in short hand notations, 

A + z; = 0, v A v = 1 N . (67) 

Note that only the lower component v 1 depends on z. 
Expressing v as 

= ( u(x) ) (p ~ m ' < X) = (B+ " X+) " lA+ (68) 

let us find u(z,x) - the main object of ADHM construction. It is the solution to the 
equation 

(B f - x f )u = A + , B + - x + = A _ r +( z ) (69) 
Define the Green's functions: 

/ = (A+A)" 1 , G = ((B-x)\B- x))~\ ^y(x) = Sij + A? G^A + J (70) 
where (i = 1 . . . N), (a, /3 = 1, 2). One can note that 

/ = (G _1 + A+A)" 1 = G - G\\(p^ l kG (71) 
acting on (71) with A + on the right yields 

GA| = fX\^ (72) 
The Green's function is expressed as follows: 

f(z,z') = S m (z)f mn sl(z') +27TS(Z,Z / )<5[ Z ][ Z '] (73) 
fmn = F mn (74) 

The functions appearing in eq.(73) are 

. (z) _ r 2mx n (z-u m ) smh[2nr m (ii m+l - z)] 2 ^ o(z _ Fm) sinhj2^^ - }i m -i)} 

sinh(27rr m v m ) ° m ^ +6 sinh(2 7 rr w _ 1 v m _i) m ^ +1 ' 



s(z,z') = e 2mx °( z -^ 



sinh(2nr [z] (mm{z,z'} - ji [z] )j sinh(27rr [z ](i/ [z]+1 - max{z,z'}) 



r\ z ] sinh 27rr fzl vr z i 



Ui = (B — x)fk}^i = (J^. - r^A s f k (z)f kj C}(t>ji (75) 

And the convenient notation is r t ■ = x — x/j is a vector from the z'-th dyon to the current 
point. First we note that Fjj and tyi are diagonal matrices with exponential precision 

fij^lnSdri + ri^ + Qi)- 1 (76) 
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* (77) 

To pass to the periodical gauge we multiply v(x) from the right by gjj = Sjje 2m l l > x °. 
Totally within the exponential precision we get for v: 

( _$.. ( p-l/2 e 2mpi i x \ 

Vi(x,z) = I , . ! V" + 1 /? t ■ no index summations (78) 

V 7 \ (B-x)s l (z)f ll Cjcpl /2 e 2m ^ J ' 

Consider the covariant derivative of Vj(x) in the periodical gauge (integration over z is 
assumed): 

D }l {v{x)\ = d }l {v\-d }l {v\v){v\ = d }l {v\{l - \v)(v\) (79) 
= d v (v\AfA f = -(rj|a F A/A + 
= -{v\Ba v f^ = (v n a }l f\\v 2 ^cr }l f(B - x) f ) 

= Hr<pl /2 e- 2m ^ X0 (B - x)'s]{z)cr }l (fA\f(B - 

= fdi<Pl/ 2 e- 2ni ^ iXo (B - x^sJizy, (sitofuG, (siMfusUzf) + 2ns(z,z')) (B - x) f ) . 

\v(x))D }l = -(D }t (v(x)) f (80) 
B.5 Formula for the regular current 

Let us denote D = (B — x). We will in a moment express eq.(59) through c and fc, 
defined as 

c = D H vlv y (81) 

= -<jmf^Q(Ds(z)Y cr }l (s(z)fs f (z') + 2ns(z,z'))e- 2nin( - z '-^Ds(z , )^e- 2m ^ n 

= -^/^■(Ds(z))V ?( (s(z)/s + (z / ) +2n~s{z,z'))D~s{z')£e- 2mnz ' 

here ~ means that the time dependence is separated (so that s is time-independent) and 
integration over z and z' is assumed. To derive this we used 

D f x v 2 = (Dj - n)^ = -Az;J - nz;2 (82) 

following from eq.(67) and noticed that the first term cancels with the scalar product 
of upper components. For b we easily get 

b = vtvy = exp(-2m> f n)/& + ^/((Ds) 1 [ e- 2 ™ z Ds)$ . (83) 
We also need the following formulas: 

v\v y D { ; j) = -cU_ n (84) 



b f (85) 



v\v x = b f 

&»— « = b f (86) 
v]v x D { y ] = -Dfv\v x + 2A }l v\v x = -c n ^- n + 2Atf (87) 

D { * ] v]v x = -v i y v x D { * ) +2Avlv x = c i + 2Ab f (88) 
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In the total sum changing n — > —n does not affect the expression since we divide by 
ft 2 , so we can make n — > —ft in the whole expression by conjugating and then drop 
n — >• —ft. 

Then the set of {(£>, C)} in eq.(59) becomes 



{(B,C)} = {(c,fc + ), (&,c + + 2A& + ), (-c^_„,& + ), (b,-c n ^- n + 2Ab f )} 



Totally we get for the current (60) (the index ji is hidden in c): 



E 

n=£0;j=\..N 



7T 2 ft 3 



/) 



(89) 



(90) 



There are only terms with (z) ^ (j) in the current that we are calculating. 
Now by these explicit formulas the current can be evaluated by performing integrals 
over z,z' and summing over ft. The reference formulas for summation can be found 
e.g. in [15]. We used Mathematica for these calculations. Below the result is presented. 
Just for illustration and consistency check we write first for the SU(2). For SU(2) the 
ADHMN data is taken to be y.\ = -co, }ii = 00 , v\ = 2co, vi = 1 — 2co. So we get for 
the current (here we write the diagonal matrix, which couples to the gauge field in the 
fundamental representation) 



r\l 



r\2 



iP' ( 27r(ft! 



. 1 1 

n) + 1 + 

r s 



12n 2 



(91) 



As usual, the current is expressed as a derivative of perturbative potential P. The sec- 
ond term in this expression cancels exactly with the contribution of "singular" current. 
Totally, adding the singular current, we get for the "far region" contribution to varia- 
tion: 



<nogDet(-D 2 ) /flr = / 5{A^ d - A^ nd )P' (ln{ H - n ) + \ - 1) 

far 

= J sp (inim - n ) + \ - 1) 

far 



(92) 



we remind the convenient notation: r l ■ — x — \j\ is a vector from the z'-th dyon to the 
current point, s, = r^_i and a standard "circular rule" r^ + i = r\. For SU(2) we set 

r — T\,S — S\ — T2- 

For the SU(N) case we get the total result: 

<nogDet(-D 2 ) /flr 



far 'i 



1 


1 






2T t 


~2s~~ 


( 2f i 


2 Sj )_ 



Sgn(ft; 



far l> i 

Note that the spatial part of the regular current cancelled exactly with the singular part 
of the current. This proves the eq.(23). 
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